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abstract 

Thermodynamics of d — 4, jV = 4 supersymmetric SU (N) Yang- Mills theory is studied with 
particular attention on perturbative expansion at weak 't Hooft coupling regime and inter- 
polation to strong coupling regime thereof. Non-ideal gas effect to free-energy is calculated 
and found that leading- and next-to-leading-order corrections contribute with relative opposite 
signs. Pade approximant method is adopted to improve fixed-order, perturbative series and is 
found to decrease free-energy monotonically as 't Hooft coupling parameter is increased. This 
may be regarded as an indication of smooth interpolation of thermodynamics between weak 
and strong 't Hooft coupling regimes, as suggested by Maldacena's AdS/CFT correspondence. 
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1 Introduction 



Maldacena's proposal [0] relates large N limit of the d — 4, jV = 4 supersymmetric SU(N) 
Yang-Mills theory to Type IIB superstring theory on AdS$ x S 5 . According to the correspon- 
dence, parameters of Yang-Mills and string theories are related each other as: 

9ym = 2 ^9st , gy M N = A. 

— LX 

where R is the curvature radii of both AdS§ and S 5 , and g st and a' are string coupling parameter 
and string slope parameter, respectively. Exploration of the proposed duality over the entire 
parameter space (gyu^ N) is seemingly a formidable, if not impossible, task. However, the 
above relation already implies that physics of super Yang-Mills theory in the strong 't Hooft 
coupling limit, A 2 = gyy^N — > oo, is directly accessible from systematic expansion in the zero- 
slope limit, (a'/R 2 ) — > 0, of Type IIB string theory, viz. gauged AdS$ supergravity. Indeed, 
all the checkpoints regarding the duality examined so far have been exclusively in this limit. 
It is therefore of interest to explore, for various physical quantities, the Maldacena's proposal 
beyond the leading order in strong 't Hooft coupling in super Yang- Mills theory and in (a' /R 2 ) in 
gauged AdS$ supergravity. Of particular interest is whether the physical quantity is an analytic 
function over the entire parameter space. If so, it would be interpolating smoothly between 
weak and strong coupling regimes. In this paper, as an example of such physical quantities, we 
would like to investigate thermodynamic free-energy density F(N, A), with particular attention 
to the issue of interpolation between weak and strong 't Hooft coupling regimes. The issue is 
particularly relevant, as, based on certain assumptions on analyticity and convergence of the 
perturbative series, it has been claimed that such a smooth interpolation is not possible ||. 

In the strong 't Hooft coupling limit, the free-energy density has been calculated in ||, 
including leading-order correction. Their calculation goes as follows. According to the Malda- 
cena's proposal extended to finite temperature, super Yang-Mills theory at high temperature 
T is described, in large N and strong 't Hooft coupling limit, by a Schwarzschild AdS$ black 
hole [|]] at Hawking temperature T. The metric of the black hole is given by 

ds 2 = H! (-f(V)dt 2 + dxj) + ^ (f-\V)dr 2 + U 2 dQl) , (1) 
where, expressed in terms of Yang-Mills paramters, the Schwarzschild factor /(U) is given by 

f(V)=U-j£j , U = 2ttAT. 

As shown in || , the metric Eq. (JJ) turns out to be a solution of gauged AdS 5 supergravity equa- 
tions of motion not only at leading-order (Einstein-Hilbert action) but also including 0(a' 3 ) 
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string correction terms. The free-energy density F is then obtained from the Euclidean super- 
gravity action of the black hole times Hawking temperature T. The result is 



F(N,T) = h(X) (-^n 2 N 2 T A 



h(X) 



(2) 



and hence indicates that the free-energy density is diminished by a factor h(X) ~ 3/4 compared 
to the free-energy of ideal gas of M = 4 super Yang-Mills theory, —tt 2 /6 ■ N 2 T 4 . Branch-cut 
of the free-energy density at A = oo is a genuine prediction for the strong coupling limit, 
along with the heavy quark potential, both at zero |5j, U and finite temperatures 0, H, ||, and 



anomalous dimensions of long supermultiplets JT0|, [I I 



Were it interpolating smoothly over the entire 't Hooft coupling parameter space, the func- 
tion h(X) in Eq.(|2|) is expected to make monotonic cross-over between h(X — > 0) ~ 1 and 
h(X — > oo) ~ 3/4. To examine whether this is the case, in this paper, we calculate the 
free-energy density of the super Yang-Mills theory in fixed-order perturbative expansion and 
examine its behavior as A is increased. According to our calculation, the leading-order 0(X 2 ) 
correction indeed makes h(X) decreasing, but the next-order, non-analytic 0(X 3 ) correction 
contributes with (relative) opposite sign. The result up to this order then seems not to exhibit 
anticipated monotonically interpolating behavior, but it may have to do with truncation of the 
perturbation expansion to lower order terms. In an attempt to improve the fixed-order pertur- 
bation theory and (partially) resum higher-order terms, we have employed Pade approximant 
method. Our result shows that, once Pade-improved, the function h(X) indeed becomes mono- 
tonically decreasing as it is extrapolated to larger values of the 't Hooft coupling parameter. 
Interestingly enough, the (relative) opposite sign of the 0(X 3 ) correction turns out precisely 
the origin of the monotonically decreasing behavior after the Pade approximant is made. 

The paper is organized as follows. In section 2, we present perturbative calculation of 
the free-energy density of d = 4, J\f = 4 super Yang-Mills theory up to 0(X 3 ). In section 
3, we improve the fixed-order perturbative result of section 2 by Pade approximant method 
and find smooth interpolation plausible. In section 4, we conclude with brief discussion on 
interpolation of free-energy density between d = 2, H = 4 supersymmetric SU(Ni) x SU(N 5 ) 
Yang-Mills theory and Type IIB supersring theory on AdS^. Preliminary results of the present 
paper has been announced by one of the authors (SJR) in ||12|| . Since then, while our paper 
is being prepared, two independent works dealing with aspects of the calculation in section 2 
up to 0(X 2 ) |13| and 0(X 3 ) [|14]] have been reported. Wherever overlapping, our results are in 
complete agreement with theirs. 
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2 Perturbative Calculation of Free-Energy 



In this section, we will be computing free-energy of four-dimensional super Yang-Mills theory 
with Af = 4 supersymmetry and gauge group G = SU(N) up to 0(<?ym) i n large- iV limit. 
Throughout the computation, we will adopt Feynman gauge and use dimensional regularization 
for ultraviolet divergences. In component fields, the J\f = 4 super Yang-Mills theory consists of 
a gauge boson, four Majorana fermions and three complex scalars in the adjoint representation 
of the gauge group G. 





/ \ 
/ \ 
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/ 
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Figure 1: One-loop Feynman diagrams contributing to free-energy F\. Gauge bosons are rep- 
resented by wiggly lines, ghosts by dotted lines, fermions by solid lines, and scalars by dashed 
lines. 



At leading order, 0(<7ym), the ^ ree energy density F\ is simply that of the non-interacting 
massless degrees of freedom and is obtained from one-loop Feynman diagrams, as shown in 
Fig.l for gauge bosons, ghosts, gaugino and three adjoint fermions, and three adjoint scalars, 
respectively. Evaluating each diagram explicitly, one obtains 

- F 1 {G,T) = -Ad G ■ ^B + 2d G ■ ^B + Ad G • 2 • ±F - 3d G ■ 2 • ^B . 

Here, d G denotes the dimension of the gauge group G, and Bq{T) and Fo{T) are defined as^ 

BAT) = %W = -^T>, 



MT - % } W = \-^. 



2 Here, [p] ({p}) represents that the four-momentum p = (po,p) is bosonic (fermionic), namely, 

d 3 p 



?fp} 



T J2 I Tfrk where Po=2nnT, 

Po 1 

r £/|^ where Po = 2n(n+ 1 -T). 
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Hence, summing up all the one-loop diagrams, one obtains Fi (G, T) as 

F^G.T) = id a (B -T ) 

D 



(3) 



The contribution of the next-leading 0(g Y u) order is calculated from two-loop Feynman 
diagrams diagrams in Fig. 0. 






(a) 



(b) 



(c) 



(d) 





(e) 



(f) 



(g) 



Figure 2: Two-loop diagrams contributing to the free-energy density F 2 



After a straightforward calculation, one finds that each diagram gives rise to the following 
contribution: 

- F 2a = -3g YM c A d G B 2 
9 

-F 2b = +^d G g YM c A d G B 2 

—F 2c = — ^9yM c Ad G B 2 
-F 2d = -16g YM c A d G {IS 2 - 2BT 

—F 2e = — —g YM c A d G B 2 

-F 2f = -\2g 2 YM c A d G B 2 
9 

-F 2g = +-g YM c A d G B 2 . (4) 
Here, c A (G) denotes the quadratic Casimir invariant of the gauge group G and 

r 1 T 2 

4 



r 1 T 2 

Adding up all the two-loop Feynman diagrams, one finds 

F 2 (G,T) = lQg 2 c A d G (B-F) 2 

= (3ym c a) tA ■ ( 5 ) 

Comparison with Eq.([3|) shows that the perturbative expansion is indeed in powers of the 't 
Hooft coupling parameter A 2 = gy M CA- Note that, at zero temperature, B = JF and B = T 
and hence bosonic and fermionic contributions in Eqs.(^), (|5|) are cancelled out, manifesting the 
supersymmetry nonrenormalization theorem for the vacuum energy. For G = SU(N), where 
da = N 2 — 1 and ca = N, the results Eqs.(|3|), (||) are in agreement with JT3], [TJ]. 

In the fixed-order perturbative expansion, the next order would be of 0{g\ M ) ) coming from 
three-loop Feynman diagrams. However, as is well-known, (a subset of) these diagrams are 
afflicted with severe infrared divergences associated with the exchange of electrostatic gauge 
bosons and scalars. Physically, these infrared divergences are screened by Debye plasma screen- 
ing. The screening of gauge bosons and scalars can be taken into account by a resummation 



of so-called ring diagrams to all orders |15|, The infinite-order resummation is the ori- 
gin of the non-analytic contribution starting at order (s'ym) 3 ^ 2 - ^ n efficient way to calculate 
the non-analytic term is to reorganize the fixed-order, perturbative expansion by rewriting the 



Lagrangian as [17 



■>SYM 



and to treat the second term as perturbative interactions. Here, M e \ is the electric thermal mass 
for the gauge boson given by the one-loop self-energy Yi^ v at zero momentum 0, M 2 l 5 ab = n^O), 
and m e \ is the thermal mass for scalars calculated from m 2 ^ * = n ab (0). Self-energy diagrams 
in Fig. 3 can be evaluated straightforwardly in the zero-momentum limit, much as was done in 
computation of Eq. The result is 

Ml = 2(g 2 YM c A )T 2 , 

m l\ = (9ym c a)T 2 . (6) 



Rather than presenting details of straightforward but tedious calculation, , we briefly discuss 

how O (54m) contribution can be obtained. The free-energy density of an ideal Bose gas with 

3 A nonzero thermal mass is generated for Aq, but not for Ai. This stems from the fact that, after dimensional 
reduction over Euclidean time direction, one still retains d = 3 gauge invariance. 
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(b) 

Figure 3: One-loop self-energy diagrams for gauge bosons and scalars. 
mass m is given by 

-Ffree boson = ~ ^/ hl(j> 2 + m 2 ) , 

= _ZL!t 4 + — m 2 T 2 -— m 3 T + C(m 4 ). (7) 
90 24 12vr V ' K J 

The first term is the well-known, black-body radiation effect of the massless scalar particles. 

Its exerts a positive radiation pressure and hence a negative contribution to the free-energy. 

The second term is the effect of nonzero mass; it contributes positively to the free energy, as 

expected on physical ground. The third term, which gives rise to a negative correction, may be 

understood as representing the bosonic nature of the particle Its contribution is of C(m 3 ), 

thus cannot be obtained in a naive perturbative expansion in powers of m? . If one interprets 

the mass m as the thermal mass due to the Debye plasma screening, Eq.(P), then Eq.(^) may 

be understood as an expansion with respect to the coupling {g\ u CA) 1 ^ 2 ■ In particular, one 

finds that the sign of the g Y M term is negative. Actually, the 0((?ym) term of the free energy 

comes entirely from Eq.(0) and there is no further correction of order g Y M- The reason is 

4 For a Fermi gas, there is no infrared divergence and hence no non-analytic contribution to the free-energy 
either. 
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simply because the thermal mass m e \ is defined in such a way that it precisely reproduces the 
resummed g\ M effect in the form of Eq.(^). Even though there are other g Y M order contributions 
from individual two-loop diagrams Fig. |2|, they are all cancelled out by the corresponding mass 
counterterm one-loop diagrams shown in Fig. 4. 




Figure 4: Mass-counterterm diagrams due to thermal Debye masses. 



Therefore, on general ground, the 0(gy M ) contribution due to Debye plasma screening 
contributes negatively to the free-energy (viz. positive contribution to the radiation pressure), 
independent of the details of the theory under consideration and can be read from Eq. (0) once 
the thermal masses are calculated. Incidentally, extra terms of order g\ M in Eq.(|7]) are also 
cancelled out by the mass counterterms given in Fig. 4, and have no net contribution, as should 
be so for consistency of the perturbative expansion. 

From the above discussion, it is now quite a simple matter to evaluate the free-energy density 
F, of 0(g§ M y 

= "^(^fV, (8) 

where the relevant degeneracy factors are also taken into account.^ Adding Eqs.@, (H) and (f|) 
altogether, one finally arrives at the free energy density of four-dimensional super Yang-Mills 
theory with Af = 4 and gauge group G up to 0(gy M ), 

F = -jd G T* [l - A ^ yuCA ) + 3 _±fi (g\ MCA f 2 + O {{g\ M c A f) } . (9) 



One can verify by an explicit calculation that the other terms of order <?ym an< ^ 5ym are an cancelled out. 
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3 Pade Approximant of Free Energy 



The next order correction to Eq.(^) will be of order g YM and can be calculated by considering 
three-loop diagrams. Moreover, a full account of the Debye screening at three loop order would 
produce terms of order ^ym^Sym anc ^ #ym- These contributions can be evaluated by utilizing 
methods developed in Refs. [|17| , [IS|. On the other hand, extending perturbative evaluation 
of the free-energy density to even higher orders breaks down beginning at order g Y M- This is 
mainly due to another source of infrared problems associated with exchange of magnetostatic 
gluons. Screening of the magnetostatic gluons cannot be achieved by resummation of the 
perturbative diagrams and, as the effect is inherently nonperturbative, can only be treated 
using a reliable method such as lattice simulation ]16) . As the new infrared divergence appears 
beginning at order g Y M5 by scaling analysis, the magnetic screening mass scale is of order gyM^- 

Rather than dwelling on estimates of these effects further, one might like to explore implica- 
tion of the result Eq.(|9|) to the behavior of the perturbative free-energy as the 't Hooft coupling 
parameter A becomes large. In particular, one would like to study if the free-energy at weak 
coupling regime Eq.@ may possibly be interpolated monotonically to that at strong coupling 
regime Eq.([|) obtained via AdS§ supergravity. On a first look, as A is increased, the positive 
sign of g Y M t erm seems to indicate that the result Eq.(|9|) at weak coupling regime behaves 
against any smooth interpolation to strong coupling regime. However, it might well have to 
do with the nature of the fixed-order perturbation theory, which truncates physical quantities 
to finite-order polynomials and induce fluctuating behavior. Certainly, if possible at all, one 
needs to improve the perturbative result Eq.(|9|) before any definite conclusion is drawn. 

Replacement of the fixed-order perturbation series by Pade approximant Jl9| may be con- 
sidered as a simple tool of gaining an idea on how the series would behave as the coupling 
parameter is increased. Indeed, the Pade approximant of QCD and other field theories, which 
has been applied already to high-energy processes |2(J quite successfully [], seems to improve 
convergence of the perturbative free-energy in the right direction We will thus apply the 
Pade's approximant to the perturbative free-energy density of the M = 4 super Yang-Mills 
theory and examine whether the approach lead to monotonic interpolation between weak and 
strong coupling regimes. 

Let us write Eq. @ as 

RW = p^t^T) = 1 + h * + /3A3 + • • • ' (10) 

where, from Eq.(|]), / 2 = — 3/27T 2 and f% = (3 + ^/2)/tt 3 , respectively. Note that, as it stands, 

R(X) exhibits a minimum at A = (—2/2/3/3), viz. Jym^ ~ 0-5- The general form of the Pade 

6 in the sense that Pade approximant " postdicts" the known higher-order perturbative terms and also reduces 
renormalization scale and scheme dependences. 
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approximants to Eq.fllOD is given by 

The coefficients c n and di are functions of fixed-order coefficients fx, f%, • • -. As the two coef- 
ficients fx and fi are determined, it is possible to construct two distinct Pade approximants, 
which reproduce Eq.fllOD up to order of A 3 : 

1 + £A 



R 



[1,2] 



1 + £A-/ 2 A 2 ' 



%i] - 1+eA , (I 2 ) 

respectively. Note that, as Eq.(^) contains no term of order A, the coefficients Cx and dx are 
fixed to the same value: 

/2 <57T 



1.5 



0.5 




I 1 1 1 1 1 

1 2 3 4 5 

(g N) 

Figure 5: Plots of fixed-order perturbative result, and Pade approximants R[i,2] and R[2,r]- 

We have plotted the two Pade approximants in Fig. 5 and have made comparison with 
the perturbative result. Note that, thanks to the fact that fz is positive, both forms of the 
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Pade approximant show monotonically decreasing behavior as the 't Hooft coupling parameter 
becomes strong. This is to be contrasted with the behavior of the fixed-order perturbative 
result, in which the strong coupling behavior is dominated by the / 3 ^y M term. Hence, within 
the simple Pade-improvement of the perturbation series, we have obtained a behavior consistent 
with the strong coupling expansion drawn from the AdS$ supergravity analysis and hence with 
smooth interpolation between weak and strong coupling regimes. 

Although the Pade approximant would not be able to draw any definite conclusion regarding 
the strong coupling behavior, at least one can say that the positive contribution of the f3g YM 
term to the free-energy density is not as disastrous as it may appear for a smooth interpolation 
over the 't Hooft coupling parameter space. Explicit calculation of higher-order g Y Mi 9ym l n 9ym 
and Qym m ight be helpful in drawing more definitive conclusion on this issue. 

4 Remark on D1+D5 and AdS^ Black Hole Free-Energy 

Finally, we would like to bring up an issue regarding another interesting configuration exhibit- 
ing AdS/CFT correspondence: a bound-state of Ni Dl-branes and N 5 D5-branes. Effective 
description of the configuration is controlled by 't Hooft coupling parameters, g Y M^i , 5ym-^5- 
At strong coupling limit, the bound-state of extremal D1-D5 brane is described by Type 
IIB string on AdS% x S3 x T 4 . An important observation of || was that the geometry is not 
modified even if the D1-D5 branes become non-extremal. This is because of the fact that the 
strong coupling expansion (viz. in inverse powers of the two 't Hooft coupling parameters) 
is expressed in power-series of Weyl curvatures and that the Weyl curvature vanishes for a 
product of AdS^ and S 3 for equal radii of curvature. This implies that the free-energy density 
of the Schwarzschild AdS^ black hole (BTZ black hole) is given, to all orders in strong coupling 
expansion, by 



viz. independent of the 't Hooft coupling parameters. 

At weak coupling limit, the D1-D5 brane configuration is described by d = 2 M = 4 
supersymmetric Yang-Mills theory with gauge group SU(Ni) x SU(N 5 ) and bi-fundamental 
matters P2| . Higgs-branch of the theory is described by four Goldstone bosons living on 
an effectively NiN 5 times bigger spatial volume. Thus, thermodynamic free-energy may be 
calculated in terms of sigma model of the four Goldstone bosons. At leading order in derivative 
expansion, the free-energy density ought to be given by 



where the ellipses denote corrections from higher-order interaction terms of the sigma model. 



F = —ttNiN^T 2 , 



(14) 




(15) 
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The leading-order term in Eq.(|15|) matches with the exact result of Eq.(14), so if the free- 
energy density were to interpolate monotonically between weak and strong coupling regimes, the 
free-energy density Eq.(|15|) should be also an exact result of perturbative expansion in powers of 
't Hooft coupling parameters. This implies that, in thermodynamics of d = 2 super Yang- Mills 
theory, there exists a curious non-renormalization theorem (not related to supersymmetry) . We 
suspect this may have to do with subtle noncommutativity between 't Hooft limit, N\,N$ — > 
oo, and thermodynamic infinite volume limit, V — > oo, in two dimensional gauge theories. 
Whether this is indeed the case would serve as an interesting checkpoint to the generality of 
the anticipated interpolation. Work along this direction is in progress and will be reported 
elsewhere. 
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